Preprint typeset in JINST style - HYPER VERSION 



Analysis and correction of the magnetic field effects 
in the Hybrid Photo-Detectors of the RICH2 Ring 
Imaging Cherenkov detector of LHCb. 



O R- Cardinale"*, C. D'Ambrosio' , R. Forty^ C. Frei^ T. Gys^ A. Petroling, D. 

^ Piedigrossi' , B. Storaci^ ' and M. Villa^ * 

^ " Dipartimento di Fisica deU'Universitd di Genova e INFN sezione di Genova, 

^"5 Genova, Italy 

^European Organization for Nuclear Research ( CERN), 
Geneva, Switzerland 

'^Dipartimento di Fisica dell'Universita di Milano Bicocca e INFN sezione di Milano Bicocca, 
''O Milano, Italia. 

[/J E-mail: [roberta . cardinaleBqe . inf n . it| 

•rH 
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^ affected by the stray magnetic field of the LHCb magnet, which degrades their imaging properties. 

O I This effect increases the error on the Cherenkov angle measurement and would reduce the particle 

identification capabilities of LHCb. A system has been developed for the RICH2 Ring Imaging 
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1. Introduction 

1.1 The LHCb experiment at the LHC 

The LHCb experiment at the Large Hadron Collider (LHC) [0] at CERN (Geneva, Switzerland) 
is designed for high precision measurements of CP violation and rare decays of heavy flavours. 
Key physics measurements to be performed with LHCb are discussed in [0]. A drawing of the 
LHCb detector is shown in Figure 0. 

Particle IDentification (PID) of hadrons is an essential tool for the LHCb physics program. 
Many final states of interesting decay channels have background from decays with the same topol- 
ogy, but different final state hadrons. 
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1.2 The RICH detectors of LHCb 

The PID capabilities of the LHCb apparatus rely on its two Ring Imaging Cherenkov detectors, 
RICHl andRICH2 OHil. 

A single RICH detector cannot satisfy the broad momentum range requirement. However 
the strong polar angle and momentum correlation can be exploited by a system of two RICHes: 
RICHl and RICH2. The two combined sub-detectors can distinguish kaons, pions and protons 
in the momentum range 1 100 GeV/c over the full angular- acceptance of LHCb by using three 
different radiator materials. 

The RICHl detector is located close to the interaction point, upstream of the LHCb magnet, 
with 25 — 300 mrad polar angular acceptance. It employs C4F10 gas and silica aerogel as Cherenkov 
radiators providing particle identification of low momentum tracks, over the range 1—40 GeV/c. 

The RICH2 detector is located approximately half-way between the LHCb magnet and the 
large ferromagnetic mass of the muon detector system at a distance along the beam axis of about 
6 m from the magnet centre and of about 10 m from the collision point. It is intended for high 
momentum, up to lOOGeV/c, and small polar angle tracks, and uses a gaseous CF4 radiator. It 
covers 10—120 mrad in the horizontal plane and 10 — 100 mrad in the vertical plane. 

Both RICH detectors have a similar optical geometry. Cherenkov photons, produced by 
charged particles traversing the radiator media above the momentum threshold for Cherenkov light 
production, are focused by spherical mirrors and reflected off flat mirrors onto a photo-sensitive 
matrix of Hybrid Photo-Detectors (HPD) \^ ^, described in Section [0[ Each RICH detector has 



two HPD planes, above and below the beam pipe for RICHl, on both sides of the beam pipe for 
RICH2. 

The RICH detectors are located in the stray field of the LHCb dipole magnet that has an 
integrated magnetic field of 4 Tm. The HPDs are affected by the stray magnetic field. As the 
HPDs must operate with optimal efficiency and imaging performance, it is important to understand, 
measure and correct for its effects [|7p. In this paper the strategy to measure and coiTcct the magnetic 
distortions in RICH2 are described. Another dedicated system has been developed for RICHl 
A drawing of RICH2 is shown in Figure 0. 




ng 



Figure 2. A drawing of RICH2 showing the beam pipe, the spherical mirrors, the flat mirrors and the HPD 
enclosure. 

1.3 The HPD of the RICH detectors 

A schematic view of an HPD is shown in Figure 0. 

Photons reaching the HPD photocathode, on the inner-side of the quartz entrance window, 
produce photoelectrons which are accelerated and focused by the internal electrostatic field inside 
the vacuum tube. Photoelectrons are detected by a silicon pixel anode chip with 32 columns x 32 
rows of pixels, bump-bonded to a binary electronics readout chip. The photocathode is biased at 
— 1 8 kV with respect to the anode. 
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Figure 3. A schematic view of the LHCb RICH HPD. 

The demagnification produced by the electrostatic field is such that the effective pixel size of 
0.5 X 0.5 mm^ on the anode chip corresponds to 2.5 x 2.5 mm^ on the HPD entrance window. 

In presence of magnetic field photoelectrons inside the HPD are deviated by the Lorentz force. 
The imaging properties of the HPD are therefore worsened. If no correction is applied, the error on 
the Cherenkov angle measurement increases, degrading the PID capability of LHCb. 

1.4 The stray magnetic field in RICH2 

The effect of the stray magnetic field on the photoelectrons trajectories depends on the intensity 
and direction of the magnetic field inside the HPD In a first approximation, when the direction 
of the magnetic field is perpendicular to the geometrical axis of the HPD, there is an image dis- 
placement; on the other hand when the magnetic field is parallel to this axis, photoelectrons move 
along a spiral trajectory causing a rotation of the image 0]. 

The stray magnetic field inside RICH2 has been calculated [jT^] and it reaches up to ~ 15 mT 
in amplitude with varying directions. In order to reduce its effects on the HPDs, the phototubes 
matrix has been installed inside a 6 cm thick iron shielding box. Moreover each single HPD is 
surrounded by a 0.8 mm thick cylindrical shield made of Mu- Metal©. This shielding effectively 
reduces the stray magnetic field inside the HPD. Simulations [jlTJ] estimate this residual field to 
be less than 1 mT with a direction approximately parallel to the HPD geometrical axis and an 
absolute value dependent on the HPD position within the shielding box. 

Measurements have been performed showing that if the field inside the HPD volume is less 
than 1 mT there is no loss of photoelectrons Qnj]. 

However this residual stray magnetic field produces a non-negligible image distortion. 

A detailed characterisation procedure of the magnetic distortion effects is necessary to imple- 
ment a correction. Measurements with magnetic field on and off have been performed, in order to 
parametrise the effects and to determine a correction procedure restoring the optimal resolution. 
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2. Experimental setup 



The correction of the magnetic distortion effect for RICH2 is based on the projection of a known 
light pattern onto the two 16x9 HPD matrices, using a commercial light projector. The pattern is 
a suitable grid of light spots and its layout can be chosen to optimise the measurements results. 

Data with the LHCb magnetic field on and off have been taken and analysed to compare the 
position of the light spots in the two conditions. 

The ultimate goal is to achieve an accuracy in the correction procedure such that the residual 
uncertainty due to magnetic distortion is negligible, in comparison to the HPD pixel size error [Q]. 

2.1 The light projector 

A commercial light projector ' [ 1I3| ] has been used to project the image onto the HPD matrix. 
Its most relevant specifications are summarised in Table |I]. 



Table 1. Main specifications of the Samsung SP-P310ME light projector. 



Resolution 


800 X 600 pixel 


Contrast 


1000 : 1 


Luminosity 


50 ANSI Lumen (1300 Lux) 


Image size at 2 m distance 


101.6 X 76.2 cm^ 


Weight 


0.7 Kg 


Dimensions (w x h x d) 


12.7 cm X 5.1 cm x 9.4 cm 



Two such light projectors have been temporarily installed during the LHCb commissioning 
phase inside the RICH2 gas vessel to shine onto the entire HPD matrix. They were centrally and 
symmetrically positioned for the two sides, at mid-height and at a distance of ~ 3.3 m from the 
HPD plane. Neutral filters have been used in order to reduce the number of photons hitting the 
HPDs, to avoid possible damage of the photocathodes. 

The choice of using a light projector, instead of a dedicated light source system, has some 
advantages. The light projector, thanks to its reduced size, is very easy to handle, transport and 
install inside RICH2. It does not contribute to the LHCb detector material budget and it is not 
necessary for it to be radiation hard. Moreover it allows to choose the shape of the projected light 
pattern with a high spatial resolution. 

2.2 The projected light pattern 

A grid of light spots on a black background has been projected onto the HPD matrix. In order 
to have the best discrimination for image movements, the smallest possible light spot size was 
chosen, corresponding to one pixel of the light projector. Since the light projector has been installed 
at 3.3 m distance from the HPD matrix, the projected size of the pixel is 2 x 2mm^ on the 
photocathode. This size is approximately equal to the RICH2 pixel size and it matches well the 
optical resolution of the light projector. 

The grid is made of single bright pixels separated by five dark pixels (see Figure This is a 
good compromise to have a large enough number of projected light spots, for better statistics, and 

'Samsung SP-P310ME 
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the need of a suitable minimum separation between two nearby light spots, to avoid superposition 
of different light spots during the image reconstruction. 
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Figure 4. One typical HPD chip (32 x 32 pixels matrix) with the projected pattern (grid of squares) and 
the reconstructed spots (superimposed dots) affected by refraction effects and reconstruction errors. Points 
on the HPD right side are not reconstructed because the HPD shielding partially shades the projector light 
which is not perpendicular with respect to the HPD plane. 

The pattern has also three larger spots corresponding to the position of three dedicated PMTs 
that have been installed in each of the HPD matrix to be used for aligning the pattern. 

2.3 The pattern alignment procedure 

Six Hamamatsu H3164-10 PMTs [jT^l have been installed inside the HPD matrix, three per side. 

These PMTs are used to align the pattern and to know the exact position of the light spots 
in the LHCb coordinate frame, since the PMT position coordinates are precisely known from the 
survey of RICH2 [Q]. This allows to check that the pattern position is always the same over time. 

In order to align the pattern, the image is moved until a maximum in the signal coming from 
the three alignment PMTs is found. In this way one can reconstruct the pattern position in the 
global coordinate system and the global position of each HPD is known. 

The data acquisition system for the alignment PMTs is described in []T3|]. 

3. Data Analysis 

An algorithm to reconstruct the position of the light spots has been implemented. It can be used 
both to reconstruct the pattern with and without magnetic field and for internal alignment purpose 
of the RICH2 sub-components. 
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Avoiding a large amount of light onto the HPD in a single shot requires to analyse accumulated 
pictures integrating on many events. The average number of hits per event has thus been set to be 
less than one photoelectron per HPD. 

3.1 Reconstruction Algorithm 

In this section the reconstruction algorithm [J7|, is briefly summarised. The algorithm is made 
of three steps: 

Search for local maxima The entire HPD chip is scanned, pixel by pixel, using a recursive 
function that looks for the absolute maximum inside a subregion suiTounding the pixel. Whenever 
a candidate maximum is found, the search region is recentered around the candidate maximum and 
the search for the maximum is repeated. If a different pixel is found with a bigger number of hits, 
it becomes the new candidate maximum, the previous one is rejected and the procedure is repeated 
again. 

The absolute maximum found in the subregion is rejected if its number of hits is smaller than the 
median value computed on the entire HPD matrix. 

All the pixels labelled as candidate maxima are then excluded and the search continues in the rest 
of the HPD chip. 

Clustering A recursive function looks for other local maxima in the region around the absolute 
maxima found at the previous step and builds the cluster. 

Estimation of the cluster center Two different methods have been implemented and compared: 
the average of the position of the pixels in the cluster weighted by the number of hits and a Gaussian 
fit of the cluster profile. They give compatible results. In both cases the standard deviation of the 
reconstructed centroids is < 0.45 pixel. 

Two main factors can deteriorate the reconstruction procedure: noisy pixels and ion-feedback 
effects: they are described below. 

3.1.1 Noisy pixel masking 

During the analysis procedure noisy pixels (< l%o of the pixels) in the HPDs have been identified 
using the reconstruction algorithm so that a masking procedure has been implemented. In fact it 
is important to apply a procedure that can distinguish between noisy pixels and genuine local light 
maxima, because a noisy pixel can satisfy the conditions to be considered as a local maximum. 
All real local light maxima, that is local maxima corresponding to a light spot, have approximately 
the same number of hits, since the projected light is uniform. On the contrary noisy pixels have a 
number of counts that is usually larger by an order of magnitude. Therefore all local maxima with a 
factor of two more counts with respect to other local maxima found in the same HPD, ai^e rejected. 
The different parameters were tuned so that noisy pixels are rejected without losing any real local 
light maximum. 

3.1.2 Ion-feedback reduction 

Residual gas inside the HPD can cause ion-feedback JTTI]; consisting in a photoelectron which 
ionises molecules of the residual gas, creating positive ions which then drift to the photocathode. 
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Ions hitting the photocathode produce secondary electron emission. A ion-feedback event is char- 
acterised by a large number of hits grouped in clusters in the HPD centre. 

For the majority of the HPDs this effect only gives a small noise increase. However for 
10% — 20% of the HPDs the ion-feedback rate can be so high that it spoils the signal and pro- 
vokes accelerated ageing of the photocathode and a decrease of the detection efficiency. 

If the ion-feedback background is large, the algorithm may have problems to reconstruct the 
light spot positions. 

A study to reduce this background has been performed analysing the behaviour of two vari- 
ables: the number of hits per event and the maximum cluster size per event. Comparing HPDs with 
large ion-feedback to HPDs without it, one can define selection cuts to reject ion-feedback events. 
A cut on the number of hits larger than three and on the cluster size larger than two has been chosen 
(see Figure The performance of the algorithm is shown in Figure ^ for one single HPD affected 
by ion-feedback. The drop of the signal efficiency due to these cuts is negligible (< 2%). 
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Figure 5. Histogram of the number of hits and maximum cluster size, used to characterise the ion-feedback 
effect in data with the light pattern (top) and dark count data (bottom). 

After removing the noisy pixels and the effects of ion-feedback it is possible to reconstruct the 
spots projected by the light projector as shown in Figure [7|. 

4. Correction of the Magnetic Distortions 

Different data sets have been taken with the LHCb magnet on and off, at full and half-field and 
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Figure 6. Hit Map for one single HPD before (top figures) and after the ion-feedback cuts (bottom figures). 
In the two top figures the algorithm is not able to reconstruct clusters due to the ion-feedback effect. After 
applying the cuts the ion-feedback is reduced and the light spots are reconstructed correctly as shown in the 
bottom figures. 



with both magnet polarities. 

A distortion effect is clearly visible in some of the HPDs. As the stray magnetic field is not 
uniform in the HPD matrix, the size of the effect depends on the position of each tube inside the 
matrix. Two examples are shown in Figure ^ 



4.1 Parametrisation 

The expected image distortion in RICH2 is a non-uniform rotation and stretch of the light spots 
with respect to the HPD centre, while translations are expected to be negligible. This means that 
the stray magnetic field has a non-negligible component along the HPD geometrical axis, while the 
component perpendicular to it is small, in agreement with simulation results jn]]. 
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Figure 7. Light spot position reconstruction procedure: raw data (top left), local reference maxima (top 
right), clusters (bottom left) and cluster centres (bottom right). 



The distortion can be parametrised by a function of the light spots position with respect to the 
HPD centre, expressed in polar coordinates {r,d): 



dtjir, 9) = AOfjir, e)ee +Ari,j{r, 0) 

where / denotes the HPD and j denotes the light spot inside it. 
The average rotation angle for each HPD is given by 



(4.1) 



(A0),- = 



(4.2) 



where is the number of light spots reconstructed in the HPD /, AOi-j is the rotation angle of the 
reconstructed light spot j with magnetic field on with respect to its position with magnetic field off 
and calculated with respect to the centre of the HPD /. 
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Figure 8. Reconstructed cluster centres for two different HPDs located in two different positions inside 
the HPD matrix : square and triangle markers for magnetic field on (opposite polarities); circle marker for 
magnetic field off. 



The average rotation angle is related to the average magnetic field {B)j at the location of the 
tube. If the magnetic field is small and uniform inside the HPD volume, one can assume that the 
rotation angle is, to a first approximation, proportional to the average magnetic field: 



(A0),oc(B).. 



(4.3) 



The sign of the rotation angle is related to the direction of the magnetic field. The effective magnetic 
field map inside the HPD volume, estimated from the average image rotation angle, is shown in 
Figure ^ The rotation angle changes its sign when going from the top to the bottom region of 
the HPD matrix and it is close to zero in the central region. The rotation effect {Ad) < 0.1 rad, is 
clearly detectable and measurable. Its value is not exactly symmetrical along the vertical direction 
since the HPD matrix is not symmetrically located inside the overall shielding box. 

A modulation as a function of the HPD columns is also clearly visible (see Figure |T^: the first 
and the last column of each side show a smaller rotation angle with respect to the central columns. 
This is consistent with the fact that the lateral columns, being closer to the vertical walls of the 
overall shielding are affected by a smaller magnetic field than the central columns. The effect of 
the overall shielding box is also evident for the HPDs installed on top or bottom of each column. 
All these effects are compatible with the results from simulations. 

Since the magnetic field intensity is different for each HPD, the distortion parameters are 
different from HPD to HPD. However, after normalising the parameters to the average rotation 
angle, which is proportional to the strength of the magnetic field, one can find universal parameters 
for all HPDs. 

Two curves are needed to parametrise the distortion effect (see Figure [TT]): the rotation angle 
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Figure 9. Effective magnetic field map inside tiie HPD enclosure, determined from the average image 
rotation angle for each HPD. Column and row numbers of the HPD position inside the matrix are given. 
The HPDs in the matrix are positioned in an hexagonal close pack arrangement. The double numbers in the 
vertical axis refer to even and odd colunms that are vertically staggered. 
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Figure 10. HPD average image rotation angle for the two HPD matrices for positive magnetic field (left) 
and negative magnetic field (right). The numbers on the axis indicate the column number of the matrix. The 
change of the magnetic field along each HPD colimm is clearly visible as well as the fact that the magnetic 
field is stronger in the central region. 
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and the radius variation of each single spot of the pattern: 

Adi-j ={Ad)ife{rrj) 
Anj = {Ad)iMrrj) 




Radius (pixel) Radius(pixel) 



Figure 11. Rotation angle (left) and radius variation (right) of the light spots as a function of the distance 
from the HPD centre. 

In Figure [TT] the quantity rj-jAdi-j was used instead of Adjj since for small r the error on the 
light spot angle is large, as distances are measured, not angles. 

From the first plot in Figure [11] one can see that the quantity rA6 increases as a function of 
the radius up to r ~ 6 pixel and then it becomes approximately constant. From the second plot in 
Figure [TT] one can assume that, on average, there is no significant change of the light spots distance 
from the centre. 

Therefore there is only a non-uniform rotation and no stretch of the image within the current 
sensitivity. 

In order to parametrise the angular distortion a fit with a second order polynomial function has 
been performed: 

fe{r)=ar + br^ (4.5) 

obtaining: 

a = 0.17±0.04pixer' = -0.009 ± 0.004 pixeP^ (4.6) 

Therefore the magnetic distortion effects can be corrected applying to all HPDs using the 
universal relations 

Aei.j={Ad)ifeirfj)/r,j ^^^^ 
An-j =0 

with the same parameters regardless of their position. Moreover the same parameters are valid, 
with a sign inversion, for the two opposite polarities of the magnetic field. 
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Measurements were also taken at half-nominal field values, demonstrating that the rotation 
angle is proportional to the magnetic field intensity. 

4.2 Correction procedure 

Once the magnetic distortion effect has been parametrised, it is possible to correct the data. The 
corrections have been applied to the light spot positions with magnetic field on and then compared 
to the positions with magnetic field off. The inverse coordinate transformations from the non- 
corrected coordinates, {x"Q^\y^Q^^), to the corrected ones, {x^'^\y^'^^) , via A0,;y, are: 

= xj"^) cos(-A0) - ') sin(-A0) 

= x^^'^ sin(-A0) + jj,'"^ cos(-A0) ^ ' ^ 

where the indices / and j have been suppressed for clarity. 

The resolution obtained after the correction procedure (see the two left plots of Figure |T^, 
amounts to 

CJv~ 0.18 pixel ay ~ 0.18 pixel (4.9) 
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Figure 12. Resolution, in the x and y coordinates, with (left) and without (right) applying the correction 
procedure (in pixels). 

This resolution must be compared to the one found without any correction of the magnetic 
field effects, (see the two right plots of Figure [T^, that is: 



RMS;, ^ RMSy ^ 0.33 pixel 



(4.10) 



For the sake of comparison the resolution due to the finite pixel size is a ~ 0.29 pixel. The correc- 
tion procedure therefore restores the optimal resolution of the Cherenkov angle [|7|]. 
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5. Correction of the magnetic distortion in the first LHCb data 

The correction parameters, determined as described before, are stored in the LHCb condition 
database and were used on the first data taken by LHCb in 2010. 

Saturated Cherenkov rings were reconstructed in minimum bias data, taken with magnetic 
field on and with both polarities. 

The rings were reconstructed with and without applying the magnetic field distortion correc- 
tion. The measured Cherenkov angle distribution of saturated Cherenkov rings was fitted with a 
Gaussian plus a second-order polynomial. The standard deviation of the Gaussian fit is summarised 
in Table 0. As the residual stray magnetic field is small in the central region of the HPD matrix, the 
data in Table show that the connection has little effect in that region. On the other hand the resid- 
ual stray magnetic field is not negligibly small in the external regions of the HPD matrix, where 
the correction improves the resolution. 

The correction parameters restore the optimal shape of the ring that otherwise would be dis- 
torted. This decreases the standard deviation of the measured Cherenkov angle distribution im- 
proving the capability of the FID system. 

Table 2. Standard deviation of the Gaussian plus second order polynomial fit of saturated Cherenkov rings 
without and with magnetic distortion corrections. 



Magnetic field down 




external 


central 


Without correction 
With correction 


0.81 mrad 
0.76 mrad 


0.73 mrad 
0.73 mrad 


Magnetic field up 




external 


central 


Without correction 
With correction 


0.76 mrad 
0.73 mrad 


0.70 mrad 
0.70 mrad 



6. Conclusions 

A system to measure the effects of the magnetic distortions in the RICH2 detector of LHCb has 
been set up. 

The residual stray magnetic field, that has a predominant component parallel to the HPD ge- 
ometrical axis, causes a small rotation of the photoelectrons from the nominal trajectory given by 
the HPD electrostatic field. 

By applying the coiTcction procedure described in this note the resolution improves from ~ 
0.33 pixel to ~ 0.18 pixel. 

The effect of the correction has been cross-checked on real data. It improves the Cherenkov 
angle resolution especially in the external region of the HPD matrix. The optimal particle identifi- 
cation capability of RICH2 is restored. 

It is planned to repeat the projector measurements during each shut-down periods to check the 
stability of the parameters and look for possible hysteresis effects. 
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